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(^ . Abstract. Lattice QCD has generated a wealth of data in hadronic physics over the 

C^ ' last two decades. Until relatively recently, most of this information has been within the 

OA , "quenched approximation" where virtual quark-anti-quark pairs are neglected. This 

f— ( ■ review presents a descriptive discussion of the effects of removing this approximation 

I— { ' in the calculation of hadronic masses. 

1 The Quenched Approximation 

In a quantum field theory involving gauge and fermion degrees of freedom, such 
^ I as QCD, we have the following path integral formalism for the expectation value 

"^ ■ of a quantity O: 

§ ■ (^> = - fvil^VipVA n{il;,iP,A) e-^^t'^-'^'^) (1) 

O' ^' 

^^ ' where Z is the usual path integral and the Euclidean action Se is defined in 

^^ I terms of the usual field strength tensor F^^, 

-)— » . 

Se^ fd''xU{x){p+m)4,{x) + ^F^,Ff,A. (2) 

^5 ■ The gauge degrees of freedom. A, are bosonic, but the fermionic degrees of 

freedom, ijj are fermionic and hence anti-commute. These are difficult to deal 
with in a computer simulation, but fortunately, since they occur as they can be 
integrated "0 analytically resulting in the usual determinant factor: 



^ 



. .C3 ; {Q) = ^ I VA Q{ii},-ip,A)det{Ip+m) e-5«('/'.^'^). (3) 

Simulations of this quantum field theory are performed on a space-time lat- 
tice by simply replacing all the continuous derivatives and integrals with finite 
differences and sums over gauge configurations. Hence, we have on the lattice: 

{n) = ^Y' Q{^,4>,A)e-^odet{4^+m) (4) 

{U} 

The (naive) lattice Euclidean action is 

S^Sf + S.^Y. {^(^+ "^)V'(a^)} - ^ E ^''('^P + ^p) (5) 

5o „ 
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where a is the lattice spacing, the Hnk variable U^{x) now carries the gauge 
degrees of freedom, and Up is the trace of the product of link variables around 
a plaquette. 

Up = Up{x)U,{x + fia)Ul{x + i'a)Ul{x). (6) 

This formalism maintains gauge invariance even on a lattice pp. 

Variations of the naive lattice action can be made to improve its convergence 
to the continuum action in two areas: 

• the naive action suffers from fermion doubling - each lattice quark flavour 
corresponds to 2"^ continuum flavours, where d is the space-time dimension; 

• lattice actions in general suffer from discretisation errors which enter when 
the continuum derivatives in (0) is replaced by the finite difference in (|5|l. 

Two methods are generally used to overcome the first difficulty ~ the Wil- 
son/clover family of actions, and the staggered action. Both of these actions can 
be tweaked so that their lattice systematic error (the second difficulty above) 
are reduced, and then they are termed "improved" . 

Simulations using the lattice formalism can be performed by replacing the 
naive sum in Q) with a Monte Carlo estimate. This introduces a statistical error 
0{1/ yjNcfg) in the estimate of (J7) where -/Vc/g is the number of configurations 
in the Monte Carlo sum. 

The lattice prescription of formulating a Quantum Field Theory has a priori 
no model assumptions - it is derived exactly from the full continuum formalism 
with no approximations. However the parameter values in real computer sim- 
ulations of lattice QCD are far from their experimental values. This is due to 
limitations in current computer power! Table Q] lists the values of the param- 
eters in typical lattice simulations along with their experimental values. Thus 
typical lattice simulations must inevitably rely on some extrapolation of lattice 
data. Note that the bare lattice gauge couphng, go, is not listed in tabled This 
because the information about g^ is contained within the a value, through di- 
mensional transmutation. Our usual intuition about high momentum transfers 
(short-distance physics) corresponding to the weakly coupled regime (small val- 
ues of go) in asymptotically free theories such as QCD, is directly applicable to 
lattice simulations. So we have go ^ as a — > 0. 

Equations Q&EJ) correctly define the full continuum theory in the limit as the 
lattice spacing a ^ 0. However, it is extremely expensive to simulate with (0]& 
0). FigureQlshows the estimated cost of lattice calculations as a function of quark 
mass using the formula in [2] for the "clover" action. (Actually, the horizontal 
axis of this plot is Mps, but, from the PCAC relation, we have Mps oc y/mq.) 
In Fig. n] we have assumed: (i) a lattice spacing of a = 0.1 /to; (ii) a lattice 
volume of L^ = (3/to)^; and (in) that there are Ncfg = 200 configurations 
in the ensemble sum in Q. (These are very conservative assumptions!) As can 
be seen, for even modest values of Mps ~ ^Mk ^ 250 MeV, full simulations 
require Tera-scale computing. ^ 

^ Recent advances in lattice actions, e.g. using an improved staggered action, have 
meant that CPU requirements are not quite so pessimistic |3- 
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Table 1. Typical parameter values in current lattice simulations 

parameter typical value experimental value 

ruq > mJlKi 50MeV mu,md^ 5MeV 

a 0.05- 0.20 /m 

L 2-4 fm oo 

Ncfg O(IOO) oo 

Nf 0, 2 or 2 + 1 "2+1" 



For this reason, H introduced the "quenched" approximation where the true 
QCD vacuum is replaced with one with no quarks present (i.e. Nf = in Table 
^. Specifically the quenched approximation is defined as follows: 

• det{^+m) is replaced by unity, thereby removing the quark-anti-quark loops 
from the vacuum configurations; 

• the coupling (3 — G/qq is shifted to try to counteract (as much as is possible) 
the removal of these q — q pairs. Typically this shift is of the following order 
pQ fa fjf'^u -(- 0.6 where /jQ'/"" refer to the coupling in the quenched and full 
theories. 

In this way the quenched approximation can be viewed as an effective field 
theory, i.e. it contains a subset of all the interactions, and the couplings of the 
quenched theory have to be tuned to take care of these missing diagrams. Figure 
12 shows the diagrams which are present in both the quenched and full theory, 
and those which are present only in the full theory. 

Quenched simulations are several orders of magnitude faster than full (un- 
quenched) simulations, and full simulations have only been performed in earnest 
in the last 5 years or so. Typical statistical and systematic errors of state-of-the- 
art full simulations are of the same order now as quenched simulations' errors 
were a decade ago. Particle physics is primarily concerned with the compari- 
son of theory with experiment, and when theoretical calculations have inherent 
errors, it is crucial to understand and quantify their scale. 

The main aim of this chapter is to determine the systematic effect introduced 
in the hadron spectrum by the quenched approximation. We will find that un- 
covering quenching effects is more difficult than one would first imagine for two 
reasons: 

• the quenched approximation proves to be surprisingly successful for many 
hadronic quantities, i.e. it reproduces much of the hadron spectrum at the 
5-10% level. Assuming that QCD is the theory of the strong interaction 
implies that removing this approximation makes a relatively small effect! 

• current full simulations have statistical errors of a few percent (since they 
are highly cpu-intensive) and so discerning the quenching effects with this 
relatively noisy data can prove difficult; 
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While the quenched approximation proves to be unintuitively accurate for 
many hadronic quantities, there are some quantities where it either fails drasti- 
cally, or has pathologies when the valence quark masses in the hadrons becomes 
vanishing. Examples of this include: 

• the J] and rj' mesons are degenerate in the quenched theory, whereas they 
not degenerate in full QCD. This is because the quenched theory excludes 
diagrams involving disconnected q — q loops. (See |H] for a description of 
lattice simulations of rj and 77'.) 

• the chiral limit of quenched QCD suffers from "chiral logs" ~ log(mq/yl;^) 
where rUq is the quark mass and A^ is a mass parameter proportional to 
pion decay constant. These logarithms enter the chiral perturbation theory 
expansion of various hadronic masses in the quenched theory, spoiling their 
chiral limit 6 

• the hypcrfinc mass splittings in heavy-mesons in the quenched theory is 
wrong by up to 10% or more, with the sign of the discrepancy depending on 
the states considered. This systematic error is greatly reduced when the full 
theory is considered. 

We will discuss some of these issues in later sections. 

The next section briefly reviews the best current results obtained from the 
quenched approximation. It outlines the accuracy of this approximation for the 
hadronic spectrum. Section 13 presents recent results from full (i.e. unqucnched) 
simulations and we attempt to uncover estimates of quenching effects in Sect.^ 

2 Results from the Quenched Approximation 

While there have been many papers published using lattice simulations in the 
quenched approximation^ we will concentrate on the work of the CP-PACS col- 
laboration who have produced one of the most accurate quenched study of the 
light hadron spectrum in 7 ^. Their calculations used an improved clover ac- 
tion PI simulated at volumes of around 2.5'^ fm^ with several quark masses and 
lattice spacings. They are thus able to perform continuum (a — > 0) and chiral 
(vTig —f rau.d) extrapolations (see Table ^. 

An impressive summary of their quenched spectrum for the light hadrons 
is shown in Fig. |21as the open symbols. Their lattice data are shown after the 
appropriate continuum and chiral extrapolations and is taken from table XV 
of [7]. As can be seen from the middle panel of Fig. O discrepancies between 
the quenched lattice results and the experimental values are 0(5 — 10%). It is 
important to note that this relatively small difference is only discernible due 

^ A search on the SPIRES database for "quenched" returns more than 500 papers, and 

this does not include papers which use quenched simulations but where the authors 

have not included this word in the paper's title! 
^ Another paper from the CP-PACS Collaboration studies even larger lattice of up 

to 64 X 112 but the lattice action employed in this work is the pure Wilson action 

which has 0(a) errors [H|. 
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to the tiny errors in the lattice data of 0(1 — 2%). Quenched calculations of 
an earlier generation |10| . with correspondingly larger errors, were not able to 
uncover deviations from experimental values. 

Figure 01 contains two sets of lattice data points: those obtained from the K 
and 4> inputs. These refer to the hadron whose mass was used to set the strange 
quark mass in the lattice calculation. (The p mass was used to define the lattice 
spacing, a.) The fact that there are differences between these two sets of data 
is itself an indication of the failure of the quenched approximation, i.e. an exact 
calculation's result wouldn't depend on how the scale was set. 

The CP-PACS collaboration also find the mass splittings, such as the hy- 
perfine splitting in the meson sector and the splittings in the decuplet (baryon) 
sector are smaller than experimental values. 

A further indication of the failure of the quenched approximation is in the 
determination of the strange quark mass m^. As mentioned above, this quantity 
can again be calculated using either the K— or (/)— meson as input, but the 
deviation between the two results is at the 3 — 4tT level. 

Moving to the heavy-hadron mass spectrum, a recent publication, using an 
improved staggered action, studied the splittings in the heavy-meson sector 
. The left-hand plot in Fig. [3 (taken from 3 ) shows the quenched predic- 
tions of various splittings from |3]. This clearly shows a discrepancy between 
the quenched results and experimental values. As we will see in Sect. 13.41 this 
discrepancy disappears when we remove the quenched approximation. 



3 Results from Full (unquenched) Simulations 

This section will give a flavour of current full QCD lattice simulations by con- 
centrating on the CP-PACS jlj and UKQCD ^J collaborations' results for the 
light hadron spectrum, and the work of [3] for the heavy-meson spectrum. Both 
the CP-PACS and UKQCD collaborations used two flavours of improved clover 
fermions whereas the collaboration in |2| used an improved staggered action 
(which has a cpu advantage over the Wilson action, see Sect. Q] and |12)). 

Table 121 displays the parameters used in these collaborations' work. Note that 
we have differentiated the sea and valence quark masses in this table (c.f. table 
^. The sea quarks are those which always appears in quark loops and are not 
connected to the source/sink operators (e.g. the quark loop in Fig. j^J, and the 
valence quarks are those which enter the source/sink interpolating operators. 

As can be seen from table (21 the CP-PACS collaboration have performed 
QCD simulations at parameter values closer to the experimental values and has 
larger statistics than the UKQCD collaboration (see also Table ^|. However the 
UKQCD collaboration chose a subset of its parameter values so that the lattice 
spacing remained flxed as the sea quark mass, m^J^"-, varied. This meant that 
0{a) effects could more readily be disentangled from dynamical quark effects. 
Furthermore, the UKQCD lattice action has the technical advantage that it has 
no 0{a) lattice systematic errors. 
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In the simulations of [3j extremely light quarks were able to be studied due 
to the use of the improved staggered action. (This seems to have become the 
action of choice for most dynamical simulations.) Also P] simulates with the 
more physical value of 2+1 quark flavours (see table QJ. 



Table 2. Indicative parameter values used in full QCD simulations by the CP-PACS 
|7| and UKQCD Collaborations |11| in their study of light hadrons, and in the study 
of the heavy-meson spectrum in ^. 



parameter 


CP-PACS IZl 


UKQCD [TT] 


Davies et al. |3| 


mf" 


0.5m,s — 1.8ma 


0.6ms - 2.0ms 


0.17ms - 0.5ms 




i.e. 0(50 - 180) MeV 


i.e. 0(60 - 200) MeV 


i.e. 0(17-50) MeV 


m^' 


0.25ms -2.1ms 


0.6ms — 2.4ms 


0.12ms — nis 




i.e. 0(25 - 210) MeV 


i.e. 0(60 - 240) MeV 


i.e. 0(12 - 100) MeV 


a 


0.09 - 0.25 fm 


~0.11fm 


0.09 fm & 0.12 fm 


L 


~2.5/m 


~1.7/m 


~ 2.5 fm 


Ncf, 


0(1000) 


0(200) 




Nf 


2 


2 


2 + 1 



Rather than give the full details of the results from these collaboration's 
work, a summary is presented in the following. In the next section we attempt 
to understand the discrepancies between this section's full QCD results and 
quenched simulations from Sect. |21 

3.1 Meson spectrum 

In Fig. 21 we plot the vector and pseudoscalar meson mass taken from [7] to- 
gether with the experimental points. In this figure, the lattice spacing, a, was 
determined from the K and K* meson masses using the method described in 
|13|. The huge number of data points corresponds to the fact that there are 
16 different {P,mg^°') combinations in [7] and that there are 9 different valence 
quark combinations for each of these {f3,m^^°') ensembles. Also plotted are the 
experimental data points corresponding to the (tt, p), (K, K*), and (77^, 0).** 

One of the main points to be taken from this plot is that the systematics 
involved in lattice simulations are clearly under control. Variations amongst this 
data in Fig. ^is ;$ 1%, while the lattice spacing and sea quark mass vary by 
around a factor of three: a ^^ 0.09 — 0.25 fm, and m**^" ~ 0.5 — 1.8™^. A close 



* Note that rjs is not a physical particle, since there is no pure s ~ s pseudoscalar 
meson due to mixing with the u, d quarks. The rja mass shown here is defined as 
M^^ = 2M]( - M^ which follows from the PC AC relationship Mpg oc m,. 
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analysis of the data has been used to extrapolate away these residual systematic 
effects in mf'^" and a [7j (see also jE]). 

Figure|2|also shows the hadronic spectrum, including the three mesons, K, K* 
and (j) taken from [Jj. These have been obtained by chirally extrapolating data 
analogous to that in Fig. 01 to the physical quark masses. As can be seen (par- 
ticularly in the lower panel of Fig.OJ the full QCD simulated results are in very 
good agreement with experiment. 

As mentioned above, the if, K* mass is used to set the lattice spacing, a, in 
Fig. 0] J3| ■ This means that the lattice data and the experimental K, K* point 
agree by construction. However, the slope of the lattice data is a real lattice 
prediction. In the next subsection we study this gradient. 



3.2 J— parameter 

In this section we analyse the gradient dMv/dMpg of the lattice data. The 
dimensionless quantity used to study this is defined ^S] 



T » ^ dMv 
dMlg 



(7) 

K,K* 



Note that we define the experimental vahie of the J— parameter by approximating 
the derivative in Q by a finite difference: 

-^ -^'^' Ml-Ml- ^^> 

Therefore the lattice estimate of J is obtained by taking the derivative in Q 
w.r.t. variations in Mpg{ing'^°',m'^"''') at fixed to^'^° but varying valence quark 
mass, 771™' (i.e. the experimental/physical sea quark mass is clearly fixed!). 

The J— parameter has been studied in both [7] and JT], but we concentrate 
here on the analysis in ^J. Figure [S] plots the J— parameter from ^Jj. The 
J— parameter is calculated at each of three m^'^" values separately. These J{m^^°') 
values are then extrapolated in m**^" oc Mpg to the physical point m^'^" ~ 0. 
This extrapolated J value is shown as a banded region in Fig. 

As can be seen from Fig. El the individual J{m'i^°') values are significantly 
smaller than the experimental value. However there is a clear trend in m^^° 
which tends towards the experimental point. 



3.3 Baryons 

The vector and pseudoscalar light-meson sector at various quark masses is de- 
fined by the plot in Fig. ^ Traditionally, the corresponding plot containing in- 
formation about the light-baryon sector (specifically the nucleon) is the "Edin- 
burgh" plot where the nucleon mass is plotted against the pion mass (with both 
masses normalised by the vector meson mass). Figure shows this plot for the 
the CP-PACS 7 and UKQCD collaborations for their unitary data, i.e. where 



8 C. AUton 

faience ^ jTj^ea ^g ^^^ t^^ seen, there is a relatively large spread in the data, 
but there is a tendency for the data to approach the experimental point as the 
quark masses decrease. 

After this chiral extrapolation is performed, the CP-PACS collaboration [Jj 
obtained the baryonic spectrum seen in Fig. |2| This is an impressive array of 
data spanning octet and decuplet sectors. As can be seen Fig. 13 the nucleon and 
A differ from experiment by around 10%, whereas the S", S* (with quark content 
Iss) and particularly the fl (with quark content sss) are in perfect agreement 
with experiment. This implies that lattice simulations become more accurate as 
the strange quark content increases [7j. 

Since lattice simulations normally have valence quarks which span the mass of 
the strange (see table I^J, the spectrum calculation of baryons containing purely 
strange valence quarks requires no valence quark chiral extrapolation. However, 
the level of chiral extrapolation required obviously becomes more and more sig- 
nificant as the light quark (i.e. u and d) content of the baryon increases. This 
suggests that chiral extrapolation procedures need more careful consideration 
in order to resolve the discrepancy above (see jE]). Note that the authors of 
Hj themselves argue that this discrepancy could be due to finite volume effects 
which impinge upon baryons composed of light quarks more than those com- 
posed of strange quarks. This could presumably also be a factor, particularly 
since finite volume effects are likely to increase the mass (which is in the direc- 
tion of the observed discrepancy in Fig.lSJ and be most relevant for the lightest 
baryons. 



3.4 Heavy Meson Mass Splittings 

There has been a recent study of the heavy-meson spectrum in [2j which uses 
2-1-1 flavours of quarks, i.e. 2 light degenerate flavours which play the role of the 
u and d quarks, and one heavier (but still dynamical) quark which plays the role 
of the s quark. This is obviously closer to the real world than the simulations of 
the CP-PACS and UKQCD coUaborations (see tables [l] & 0) . 

We reproduce, in Fig. d (taken from 3;) a graph showing the ratio of lattice 
prediction to experiment for some heavy-meson mass splittings. As can be seen, 
the lattice results in the full theory (right-hand plot) are within la of their 
experimental values. 



4 Quantifying Quenching effects 

4.1 Hadron Spectrum 

In Sects. 13 and 13 we have outlined some results for the hadronic spectrum for 
both the quenched approximation and full QCD. Comparing these results we 
note firstly that quenched results are generally within 10% of their experimental 
value for a wide variety of quantities. This is an unexpectedly good level of 
agreement which will be discussed later in this section. 
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Studying the light-meson sector, we note that the fuU theory is able to ac- 
curately reproduce the K, K* and (j) masses to within around 1%, far better 
than the quenched theory (see lower panel of Fig. |2I). We note, however, that 
the "chiral" slope, defined via the J— parameter, (0, is still several a away from 
its experimental value at the simulated values oirrisea (see Fig.EJ, and that the 
chiral extrapolation, nisea —^ 'mu,d is required to make contact with experiment. 

This situation is mirrored in the baryonic spectrum. Figure |21 shows the re- 
markable prediction from the CP-PACS collaboration of eight baryonic masses. 
In general terms, the agreement between theory and experiment is enhanced 
when the quenched approximation is removed. Note also that there is no dis- 
crepancy in the full theory between predictions using the K and mesons to 
set the strange quark mass. The same is not true in the quenched data (see Fig. 
01 . The level of agreement between the full theory prediction and experiment 
is most profound for baryons containing the largest strange quark content. We 
argued in Sect. 13.31 that this could imply that the lattice data at the simulated 
values of m^ (roughly around nis) are correct, but that the chiral extrapolation 
procedure iriq -^ rriu^d is going astray. As can be seen in Fig.|Hl which is roughly 
the baryonic equivalent of Fig. 01 the chiral extrapolation required to reach the 
u, d quarks is substantial. 

Moving to the heavy-meson sector, we summarised in Sect. |^3 results from 
0. These show excellent agreement between full simulation results and experi- 
ment for a variety of quantities, especially splittings in the v spectrum. A cor- 
responding quenched analysis shows discrepancies of ~ 10%. 



4.2 Why is the Quenched Approximation so good? 

While there are obvious failures in the quenched approximation's ability to re- 
produce the real world, it does much better than naive expectations: one would 
imagine that removing all q — g diagrams from the theory would have a drastic 
effect on the hadron spectrum. Figure O shows that this is not the case. Why 
then does the quenched approximation perform so well? 

One can obtain a handle on this issue by studying the static quark poten- 
tial (which is the quantum mechanical potential between two infinitely heavy 
quarks). Figure |H1 shows UKQCD results for this quantity for both the quenched 
and full theory 11 . The curve shown in the graph is the "string model" , V{r) = 
e/r + a jr + const | 17| . Note that the data is defined to agree in value and slope 
exactly at r = tq (the hadronic scale defined in TS") TT". A close up of the 
difference between the lattice potential and the string model at short distances 
is shown in Fig.|51 As can be seen from Figs.|Hl&ini the discrepancy between the 
quenched and full theories is negligible across the whole range of r except at very 
small distances where the deviation is discernible, but small. This implies that 
only physical quantities particularly sensitive to this short-distance scale will be 
affected by the quenched approximation. Hadronic states are most sensitive to 
"medium" distance scales r ^ rg Ki 0.5 /to where (from Fig.lHJ the two theories' 
data overlay each other. 
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Thus the quenched and fuU theories should agree at the same level as quarks 
in QCD can be approximated as moving in a static quark potential. This obser- 
vation presumably has relevance to the age-old question: Why does the (non- 
relativistic) quark model perform so well? 

It is worth noting that, from Sect.^ the quenched approximation is defined 
not just by replacing the quark determinant by unity, but also by renormalising 
the coupling go. In fact, if you attempt to perform quenched and full simulations 
at the same value of go, then the lattice spacing, a (or equivalently the cut-off 
'^ 1/a) will differ by a factor of around four. This is telling us that the virtual 
quark loops really are affecting the dynamics of the simulation. The apparent 
contradiction between this fact, and what we have seen above, i.e. that the 
quenched approximation reproduces the full theory (at the 5 — 10% level) is 
resolved as follows. The lattice only actually predicts dimensionless quantities, 
normally expressed as, e.g. M x a, where M is some mass. In this way the lattice 
is able to predict dimensionless ratios of physical quantities only, e.g. M1/M2. 
Although switching the quark determinant on and off does directly affect the 
lattice spacing, and therefore Ma, it seems to have little effect on the ratio 
M1/M2. In other words the physical prediction from the lattice of AIi, which 
can be obtained from M1/M2 x Mg^^ , doesn't seem to be greatly affected by 
quenching. This is telling us something remarkable: for a wide variety of hadronic 
masses (and the static quark potential), the removal of virtual quark loops from 
the theory can be counter-balanced simply by an adjustment in the coupling, 

go- 

There is one final reason why quenching has only a modest affect on the 
hadronic spectrum, compared to other physical quantities. In order to extract a 
hadron mass from a lattice simulation, the quantity J7 = C{t) is calculated (see 
(0J) where C{t) is a two-point correlation function between hadronic currents. 
In Euclidean space-time, we have 

C{t) ^ M'^e-^'^* as t ^ 00 (9) 

where Al is a matrix element between the vacuum and the hadronic ground 
state. Lets assume that we are performing a quenched calculation of C{t) and 
that it has a relative error of e due to this quenched approximation, i.e. 

C(t)'5 =C(i)^""(l-Ke), (10) 

where Cit)^'^"^^^ are the quenched and full correlation functions respectively. 
Because the mass, M , appears in the argument of the exponential, a relatively 
small adjustment in M can mop up the quenching error e, whereas a larger 
relative change would required of the matrix element, AA. 

Obviously this analysis is a little simplistic but it does illustrate that we 
can expect quenching errors in matrix elements (such as decay constants) to be 
larger than in masses. 
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5 Conclusions 

Lattice QCD is an approach to solving field theories, such as QCD, which involves 
no model assumptions. Given a fast enough computer, the lattice can be used 
to solve QCD on any finite volume and with any non-zero quark mass resulting 
in an absolute theoretical prediction of QCD. However, in order to make the 
problem tractable on current computers, certain parameters of QCD need to take 
non-physical values (see Table ^ . The parameter under study in this chapter 
is the number of quark flavours in the vacuum, Nf. Setting Nf = is called 
the quenched approximation and corresponds to ignoring virtual q — q pairs 
in the vacuum. The approximation Nf = is seemingly a particularly brutal 
approximation and, furthermore, there is little theoretical guidance as to its 
effect. Thus we are usually forced to "measure" its effect a posteriori by analysing 
data from lattice simulations. 

In this chapter we have studied the hadronic spectrum with and without the 
quenched approximation, in particular light mesons and baryons, and heavy- 
mesons. By comparing quenched data with experimental masses, we have shown 
that quenching effects in the light-hadron spectrum are relatively small (5 — 
10%), with a slightly larger discrepancy in the heavy-meson spectrum. It is 
only recently that full QCD simulations have been able to produce data with 
statistical and systematic errors beneath this level. With this new generation of 
data, we are now able to state that full QCD lattice results have better agreement 
with experiment than quenched results. 

We have outlined some reasons why the quenched approximation is so rel- 
atively successful, and we have found evidence that the chiral extrapolation 
techniques currently being used in full QCD simulations require further consid- 
eration. 

In the future, more precise calculations with, and without the quenched ap- 
proximation will surely enhance our understanding of the underlying physics of 
QCD. 
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Mp3 [MeV] 



Fig. 1. The computer time in Teraflop-years required for a full lattice QCD simulation 
as a function of pseudoscalar meson mass using the formula for clover actions in |5|. 
We have assumed (i) a lattice spacing of a = 0.1 fm; (ii) a lattice volume of {3fm)*; 
and (iii) that 200 configurations in the ensemble sum in JIJ are required. The physical 
points corresponding to the vr— and K-mesons are shown by vertical lines. 
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a) 



b) 



Fig. 2. Diagrams wliicli are present in (a) botli quenclied QCD and full QCD, and (b) 
present only in full QCD. The full lines are quarks and the spirals are gluons. 
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Fig. 3. The Hght hadron spectrum Irom CP-PACS [3, tables XV and XII. Both the 
quenched and lull QCD results are shown, together with the experimental value. The 
two hadrons on the left do not contain strange quarks, whereas the other hadrons 
do. The lattice spacing was set from the p— mass. Two methods were used to set the 
strange quark mass: from (i) the K and (ii) the 0— mass. The results of both definitions 
are shown. In the middle and lower plots, the relative deviation {— Miatuce/Mexpt) is 
shown. The lower plot is a close-up of the middle plot showing the relative deviation 
for the strange mesons. 
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Fig. 4. The light meson masses from CP-PACS (see also IT^ '). In this figure, we 
have set the lattice spacing from the method described in |13| . Also shown are the 
experimental points. 
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Fig. 5. J versus (Mp^''"'^'')^ using the approaches described in the text from 
|11|. ]^^*^°-^y js ^}je pseudoscalar meson mass comprising of degenerate valence 
quarks which are themselves degenerate with the sea quarks. (Note, from PCAC, 
{M^ ari;-j2 ^ rUq'^"' .) The quenched data points have been plotted at (A/pg* aryy. _ 
1.3 GeV^ for convenience. The banded region at the left of the graph is the result of the 
extrapolation m''q°' -^ for the full QCD data. The experimental value of J = 0.48(2) 
is also shown. 
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Fig. 6. The "Edinburgh plot" for selected full QCD data from |7ITT] . The lattice data 
points shown are the unitary points (i.e. tti^"''^"'^'^ = rriq''"'). The experimental point is 
shown, along with the static limit {niq — > oo). 
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Fig. 7. Heavy-meson mass splitting (together with some light hadron quantities) taken 
from 13 • 



Unquenching Effects in Lattice Hadron Physics 



19 



2.5 



1.5 



> 

S -0.5 



-1.5 



-2.5 



M' 



,<v. 



• 5.93 Quenched, 


623 


- 5.29, c=1 .92, k= 


.13400 


□ 5.26, c=1.95, k= 


.13450 


<> 5.20, c=2.02, k= 


.13500 


' 5.20, c=2.02, k= 


.13550 


+ 5.20, c=2.02, k= 


.13565 


- Model 





r/rO 



Fig. 8. Tlie static quark potential Irom the UKQCD Collaboration |11|. The param- 
eters c and k refer to a coefficient of an improvment term in the action and the sea 
quark mass parameter respectively. 
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Fig. 9. The deviation of tlie static quarlc potential in Fig. |H]from the string model 
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